Complex Analysis of a Piece of Toda Lattice^ 



Satoru SAITO, Noriko SAITOHt, Hisao KONUMA and Katsuhiko YOSHIDA 

Depertment of Physics, Tokyo Metropolitan University, 
Hachiohji, Tokyo 192-03 Japan 

^ Department of Applied Mathematics, Yokohama National University, 
Hodogaya-ku, Yokohama 240 Japan 



Abstract 
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1 Introduction 



The two dimensional Toda lattice is one of soliton equations which has become more and more important 
as a key object in theoretical physics. It was first formulated by Hirota in 1981 [|] as a discrete version 
of two continuous time Toda lattice and shown by Miwa|^ its equivalence to the KP hierarchy. When 
quasi periodic solutions are substituted, it is nothing but the identity known as Fay's trisecant formula 
which characterizes algebraic curves. 

This equation has become known in other fields of physics in the last ten years. It was shown 
being satisfied by the string amplitudes in particle physics]^. More recently there appeared papers 
demonstrating unexpected correlation of this equation with other topics in physics. The transfer matrix 
of the solvable lattice model with Ai symmetry, for example, was shown to satisfy this equationQ]|^. 
This equation has been also proven to unify discrete Painleve equations The connection of solvable 
cellular automata to this equation offers another example H. 




Completely integrable nonlinear systems must play fundamental roles in various phenomena in physics. 
It is remarkable that many integrable systems in different fields are unified into the single equation. We 
are interested in clarifying ultimate notion of intcgrability of the systems. Investigation of such systems 
themselves, however, will not reveal all of features of the systems. The real meaning of intcgrability will 
be clarified only in comparison with nonintegrable systems. 

An arbitrary deformation of the two dimensional Toda lattice will destroy intcgrability and create 
chaos. Since the system contains infinite number of degrees of freedom it is extremely difficult to study 
analytically the behaviour of transition from nonintegrable to integrable phases. It should be recalled 
that a very little is known about analytical properties of nonintegrable systems. The main part of the 
studies of complex dynamical systems were limited to simple systems with one degree of freedom. 

Very recently we pointed out|^ that a set of lattice points, which form a parallelogram in the two 
dimensional lattice space, constitute a piece of the Toda lattice. We call it a Toda molecule^ since 
it is essentially what is intended to be called by this name, but used in a bit different context in the 
literature. The remarkable fact is that the small pieces can be separated from other parts without loosing 
any properties of the original Toda lattice. 

The purpose of this paper is to study in detail analytical properties of the smallest piece of Toda 
molecules. The smallest Toda molecule is a smallest parallelogram of four lattice points. We will call it 
a Toda atom for convenience. Since every Toda molecule preserves properties possessed by the Toda 
lattice, we can study analytical properties of the system from the knowledge of a Toda atom. In the first 
part of this paper we show that the time evolution of a Toda atom is described by an iterative Mobius 
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map. The form invariance of this map certificates integrabihty of this system. In the second part of this 
paper we will consider a deformation of this piece. Under generic deformation a chaos will be generated 
through the time evolution. We are especially concerned with analytical property of the Julia set as the 
system approaches to the integrable map. We will show how the Julia set converges to the points on the 
orbit of Mobius map as a parameter, which interpolates between integrable and nonintegra 

2 Pieces of Toda Lattice 

In this section we like to show that the two dimensional Toda lattice can be cut into small pieces without 
loosing any properties possessed by the original system. To begin with let us write down the equation 
which was derived by Hirota as a discrete version of the two continuous time Toda lattice 

a gn{l + l,m)gn{l,m + l) + f3 gn{l,m)gn{l + l,m + 1) 

— (a + /?) g„+i(^ + 1, m)g„_i(Z, TO + 1) = 0, a,/3eC, l,m,neZ. (1) 

We called this equation Hirota bilinear difference equation]^ and abbreviated as HBDE. This is a nonlinear 
system defined on the three dimensional lattice space. Our key observation is the following. For a fixed 
point of the lattice (/, m, n) — {I, rh, n) , we denote by A the set of points (I, fh, n), (l + 1, to, n), (Z, to + 
1, n), (l + 1, TO + 1, n), (F + 1, TO, n + 1), (r, TO + 1, n — 1) . Then if g,i{l, to) is a solution of (|^), 

f(l ^ ^ f 9n{l,m), {l,m,n)eA, , . 

/(i,TO,nj-|p^ otherwise, ^''^ 

is also a solution of (|l|) . This is the smallest piece of the Toda lattice. 

The proof is simple. Because A is surrounded by zero, every equation on other pieces is automatically 
satisfied. The result can be easily generalized to larger parallelogram prism when it is surrounded by 
zero. We call it a Toda molecule according to ref. Q. Then it will be natural to call (|^) a Toda atom. If 
there are many Toda molecules in the three dimensional lattice space separated by zeros from each other 
it is again a solution of (^. A slice perpendicular to the I axis of such example is presented in Fig. 1. 

^This equation is also called Hirota-Miwa equation in recent literature. 
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For an illustration let us consider the one soliton state localized on the smallest parallelogram specified 
by (to, n) — (0, 0), (0, 1), (1, —1), (1, 0) on the (to, n) lattice plane, but allowed to range all integers along 
I. Now we recall that in the usual lattice space the one soliton solution is given byfl [p^ 

m, n) - - az,)-''^ + I](l - bz,)-"^' . (3) 

j 3 

Here o, b are arbitrary constants and {zj\ are parameters which determine velocity of the soliton. {fc^} 
are variables taking values on integers. We can choose any three among {%} to relate them to our 
variables (Z, to, n). Let ki,k2, fcs be such three and relate them according to 



fci=TO + n-i, /c2 = -?7i-^, fc3 = Z-n-i. 



(4) 



Writing (|^) explicitly we find 

f '"'{1,0,0) 

/!-'(?, 1,0) 

f'°'il,0,l) 

/!-'(?, 1,-1) 

where 
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— azi 
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- az2 



(1 - azs)-' + B 
{l-az3)-'+B 







1 - 


bz2 


1 - 


bzi 


1 - 


bz-i 


1 - 


bzi 


1 - 


bz2 



(1 - hz^r' 
{i-bz^r' 



1 — az3 



1-6^3 



(5) 



A v/(l-azi)(l-az2)(l-az3), B ^J {I ^ bzi){l ^ bz2){l - bz^). 

We see from (||) that all points belonging to the same piece behave similarly. The parameters are related 
to a, (3 of (pi) by 



a = zi{z2 - Z'i), 13 = Z2{z3 ~ zi), 



for (ph to satisfy HBDE. If we define the amplitude (fimnil) by 



f{l + i,m,n)f{l-l,m,n) 



it behaves as 



cosh2(p? + x)' 2 I-6Z3' 2 



This represents a locahzed peak along the / axis. The other amplitudes </3mn(0 behave almost the same 
but different by the values of the phase x- 

If we consider an evolution of the system in variable /, a Toda atom is composed of four lattice points. 
Since there is only one equation of motion they are not independent variables. Three of them can 
be chosen as we like leaving one to be determined by the equation. Let zi be /(/,0,0). The other three 
could be either dependent or independent of zj. If they are independent of zi, the equation of motion 
is linear in zj. On the other hand if they do depend on zi they are allowed at most linear in z;, for the 
equation to remain Hirota bilinear form. Namely we can write 

f{l,m,n) = A„,^r^Zl + B,n,n, (m, n) = (0, 0) , (1 , 0) , (0, 1 ) , (1 , - 1) , (8) 

with ^0,0 = 1, Bo,o = 0. Upon substituting them together into HBDE, it is easy to see that we obtain 
an equation of the form 

Azi + B 

A = -/3Bi,o + {a + /3)Bo,iAi,_i, B^{a + /3)Bo,i5i,-i, 
C = {a + /5)(Ai,o - Ao,i^i,-i)- D = aBi^o - {a + /3)Ao,i^i,-i- 
is a Mobius map. Thererfore the map is integrable. 

If we remember that HBDE is invariant under the transformation of f{l, m, n) ^ai+bm+cn ^^-j^ 
the one soliton solution (||) offers an example of (||). 

The solution of can be obtained as follows. A Mobius map has three fixed points. By an appro- 
priate transformation : — > o z; o c/;^^, one of the fixed points can be transformed into 0. After the 
transformation the map will have the form 

zi+i^H—^ — . (10) 
( |l0|) is easily solved for an arbitrary initial value zq to get 

1 + ^T^^O 
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Applying to this the inverse transformation : — > o z; o , the general solution to (|^) is obtained. 
We call the map ( pX{ ) the integrable logistic map (ILM). The meaning of this name will become clear 
later. 

We notice that (^o|) corresponds to the case in which one of the lattice point is fixed constant and the 
other three points behave the same: 

f"^'' (1,1,0) = M = (12) 



How does the amplitude look like in this case? To see it we substitute ( pT[ ) into and get 

ILM sinh^p 1 1 ^zo ,.„x 
= , 2. , , — ^ n-' P:=oln^, x:=oln:j — — — • (13) 



The similarity of this result to the one soliton solution must be apparent. 
We may further simplify the equation by 

/™(/,0,0)=:zi, /'-(/,l,-l) = /™(/,l,0)-l-i, /'™(Z,0,l) = c(c: const). (14) 

A* 



The map turns to be linear 



and yields the solution 



The corresponding amplitude is 



zi+i = fizi + (1 - fJ.)c (15) 



=^'(zo-c)+C. (16) 



iinn\ sinh^p 1 1 Zo - c 

^ (^^ = cosh-(,^ + x) ' ^=2^^'^' ^=2'-^- 

which is again the form of (^. 

3 Generalized Logistic Map 

As we have learned in the preceeding section the smallest piece of Toda lattice already possesses useful 
informations of the integrable dynamical systems. In this section we study a deformation of the Toda 
atom. There could be many different ways of deformation, some of which preserve integrability and some 
others destroy it. Since we are interested in studying the transition between integrable and nonintegrable 
maps, we must break integrability of the Toda atom. 

For this purpose we recall that the Toda molecules have a characteristic form as seen in Fig. 1. 
Their cross sections in the (m, n) plane are parallelogram declined to the same direction. It owes to the 



property of the Toda atom defined in (j^) . The very reason of this asymmetry comes from the asymmetry^ 
of HBDE under the exchange of I and m as seen in . The equation in which the role of I and m in (|^) 
are exchanged is also integrable. In fact we could start from it without changing none of the results. 
From this argument we are tempted to consider the following deformation of HBDE. 

a g„{l + 1, m)gn{l,m+ 1) + /? g„(/, m).g„(Z + 1, m + 1) 
- (a + /3) [(1 - j)Sgr,+i{l + 1, m) + 7g„+i(/, m + 1)] 

X [{l-j')5'gn-i{l,m+l)+-f'g„^i{l + l,m)]=0. (18) 

We notice that this equation is integrable when 7 = 7' = and (5(5' = 1 ,or 7 = 7' = 1. Integrability of 
other cases is not known at this point. Moreover we are not able to separate some small part of lattice 
independently from the rest as it was done to get a Toda atom. Nevertheless it is worthwhile to study 
(|l8|) defined on a portion of the lattice shown in Fig. 2a. 

Fig. 2a Fig. 2b 
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In order to proceed further we have to specify the model so that we can study analytical properties 
of the map explicitly. We will consider, in the following discussion, the map given by ( |ic| ) and its 
deformation. We also restrict our argument to the case of 7' = 0, (5 = 5'^^ = ^ in ([l^ ) for simplicity and 
define (Fig. 2b) 

f^^'^il, 0, 0) = f^^'Hl, 0, 1) = /^^*'(;, 1, -1) = /^^*'(/, 1, 1) =: 

/°^*^(Z, 1, 0) = '-^-^ {^l : const). (19) 

V 

The dynamics of this model is described by the map 

l + v[l- j)zi 

We call this map a generalized logistic map (GLM). Some properties are listed below : 

■^If we had chosen other set of variables, HBDE looked more symmetric and the corresponding Toda 
atom could be either cubic or octahedron[|ll|. We have used asymmetric variables such that deformations 
can be discussed. 
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1. When 7=1, the map becomes the ordinary logistic map studied in the literature extensively. 

2. GLM becomes the logistic equation for all values of the parameters 7,/i, when the continuous 
limit of the variable I is taken. To show it let us introduce new variables u and new parameters a 
and h by 

u{l) :— zi, an := ji — 1. (21) 

^i- 1 

We replace zj+i by zij^h and take the limit ft, ^ 0. We will find that ( po[ ) reduces to the logistic 
equation: 

du , ^ , ^ 

— =au(l-it). (22) 

3. GLM includes ( p^ ) as the special case with 7 = 0. This explaines the name of ILM used for (|l^). 

4. GLM generates Julia set as long as 7 7^ 0. Hence it is not integrable except for 7 = 0. This will be 
discussed later. 

The most important feature of GLM is that it interpolates nonintegrable map to integrable map in 
the limit of continuous deformation. This fact enables us to study analytically the transition between 
two phases. The problem we concern in what follows is the analytical properties of the map (po|). To 
proceed further it is more convenient to convert the map ( ^o|) into the standard form of rational map of 
degree 2: 

F(z) = 00/0 (z) = iili-^ e»« (23) 

1 -\- A Z 



by the Mobius transformation 



where 



'^^ ^ {v^ — v — ^'^)x + {v^ — V — ^)lj,e ^ ^ 



X = ^e-^'>, ^Y2_J^^ln^,^r (25) 



The corresponding integrable map turns to be the following case 

F{z) = fiz = Xe'^'z, if 7 = i.e. AA' = 1. (26) 

The main feature of a dynamical system is determined by the nature of fixed points of the map. 
Namely the multiplier A at fixed point a of the map ip{z) is defined by the derivative of the map at a: 

d(p(z 



A := 



dz 

The fixed point a is an attractor of the map if |A| < 1, a repeller if |A| > 1, and neutral if |A| = 1. 



(27) 
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In the case of GLM, the fixed points are easily found as 

1 — AA' 

0, p=-X-^^^—^, oo. (28) 

The corresponding multiphers are 

Ao = Ae^^ Ap= ' l^y^^, , A^^A'e-^ (29) 
In the integrablc hmit AA' ^ I, we observe the following characteristic features: 

1. Since 

lAoA^I ^ 1, (30) 
the map converges either to or to oo depending on |A| — \fi\ < 1 or > 1. 

2. The fixed point p approaches to — A and it turns to a super repeller 

|Ap| ^ oo. (31) 

4 Julia Sets 

In the complex dynamical systems, chaos appears from a Julia set. Given the map f{z) on a Riemann 
sphere C = C U {oo}, the Riemann sphere is devided into two parts depending on whether the orbits 
converge or not. A set of initial values whose orbits, together with their neiborhood, converge is called 
Fatou set F{f). On the other hand, a set which does not is called Julia set J(/). This definition leads to 
the fact that the Julia set does not contain any attractive periodic cycle. In this sense the orbit in Julia 
set is chaotic. 

By definiton, Fatou set and Julia set are invariant of the map, that is 

f{F) = r\F) = F, /(J) = r\j) = J- 

It is easy to understand that attractive fixed points belong to Fatou set. Contrary it is known that 
repulsive fixed points belong to Julia set jl^ . Then we can compute Julia set by inversely mapping a 
repulsive fixed point as an initial value. We show some of their examples in Fig. 3 for the map of (p3|). 

The Julia set does not exist if the map is completely integrable. Integrable maps converge to orbits 
predictable for any given initial values. Conversely if there exists an orbit not predictable for some 
initial values, the map is not integrable. Therefore a Julia set appears in nonintegrable maps, but not in 
integrable maps. 

In our standard map of degree 2 given by (^3|), a Julia set is known to exists except for at the 
integrable point AA' = 1 . We like to know how it disappears from the complex plane of the variable when 
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the parameters approach to the hmit AA' — > 1. We have given in an argument about this problem 
for some limited range of parameters. The purpose of this section is to present another argument which 
should supplement our previous one. 

The inverse map of (^3|) is easily obtained as 

zi - F-i(z,+i) = - A) ± \sj{pzi^x + A)2 + 4zi+ie-*»(l - AA'), (32) 

where we defined 

p ■= A'e-*^. (33) 

From this expression it is apparent that the inverse map is not unique but double valued at every step. 
As we pointed out in above the inverse map generates points of the Julia set if it starts from a point on 
the Julia set. Substituting one value of the Julia set into (|3^), we get two points every time. After n 
steps the number of points of the Julia set increases as many as 2""''^ — 1. This explains the nature of 
the Julia set. Some of the points could be those of periodic maps. They must be subtracted from the 
number. 

In the integrable limit AA' — > 1 the inverse map (|3^) is still double valued. They are 

We notice that the second solution does not depend on zj+i, hence is the same at every step of the map. 
For ( p^ to generate the Julia set we must start from a repulsive fixed point. When |A| > 1, the origin is 
such a point. Thence we find from (jsj) the 'Julia set'^ 

J^^^^ = {-p"A I neiV} (35) 

for the integrable map. The number of the 'Julia set' increases proportional to the number of the steps 
n. Moreover the element of J^^^^ is equal to —fj^^^X, which is nothing but the solution exactly expected 
from the map (p6|), if it started from —A. 

The next problem we concern is to explore how the Julia set of GLM turns into those points of (jsj) in 
the limit AA' 1. Since we are interested in the transition from a nonintegrable map to the integrable 
map, we are to consider small values of |AA' — 1|. The inverse map ( |3^ can be rewritten as 

F-i(z) = |^^^|±i?(z) (36) 

where we put 



^This set does not possess properties expected for the ordinary Julia set. We call it 'Julia set' only in 
the sense that it is generated by the inverse map starting from a repeller. 
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Note that E{z) vanishes for smaU values of e. To see the behaviour of E{z) for small e we first observe 
the inequality which is true for all e: 



\E{z)\<3y/\Ie\, ^eeC. 



The proof of this inequality owes to the following facts. 



1. If \w\ < 1, 



(Vl-w-l) 



1 



1 



|l - Vl - w\ < 



< (i-(i-|^|)) = |v^|<i. 

Iv^l 



2. If \w\ > 1, 



'w ^ ' 



w \ w 



< 3. 



Substituting 



into E{z) of (|37|), we can write 



Azee 



(pz + A)2' 



\E(z)\ = Jlzel 



(Vl - w - 1) 



from which (^8|) follows. 

We can perform the inverse map (|36|) iteratively. Let us denote the map (p 



as 



A{z) := pz + E{z), B{z) := -X- E{z). 



Then the second map becomes 



F-\z) = 



A{F-\z)) ^ 
B{F-\z)) 



AoA{z) 
AoB{z) 
B o A{z) 
B o B{z) 



After n steps we obtain 



F-''{z) = {A''^oB''^oA'- 



i?"" (z) I Vi + V2 + ■ ■ ■ + Vn — n} . 



(38) 



(39) 



(40) 



(41) 



(42) 



(43) 



(44) 



(45) 



If we had started from the repeller the above maps have produced the Julia set of GLM. In the 
following we consider the case |A| > 1, |A'| < 1, so that the origin is a repulsive fixed point and the 
infinity is an attractive fixed point. Since E{Q) — the origin is mapped to 



A(0) = 0, B(0) = -A 



(46) 
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by the first iteration. The second iteration yields 

^2(0) = 0, ^ o B{0) = -p\ + E{-X), B o A{Q) = -A, ^^(0) = -A - £;(-A). (47) 

We notice that since E{—X) is the order of e from (^^, all of the points after the second iteration are in 
the neighbourhood of J^^*^. Proceeding similarly we obtain the Julia set as follows: 

jGLM = { A"' o B"^ oA""" 0---0 B"^ (0)1 1/1,^/2, • • • e AT} . (48) 

We remark some important properties which result from this expression. 

1. The invariance of the Julia set under the map. 
It is obvious from ( ^8| ) that 

jGiM = A (JGLM^ u B ( J'^^^) = ( J«^*^) . (49) 

2. An element of the form B o X for any X e J^^^^ belongs to the neibourhood of —A, as seen from 

BoX ^-\-E{X). (50) 

3. An element of the form A'^ o B o X maps _B o X to the neighbourhood of —p^X. 
In fact after applying A's s times we get 

A" (BoX) = A'-\pB oX + E{B o X)) 

= A''-'^{p'^BoX + pE{BoX) + E{AoBoX)) 

= p''BoX + ^p^E{A''-^-^ oBoX) 

k=0 

s-1 

= -p'X- p'E{X) + ^p^E{A''-^-^ oBoX) . (51) 

fc=0 

Since every element of J*^^^^, beside 0, is either the form oi Bo X or A'' o B o X, we conclude that every 
element of j'^^^^ is in the neighbourhood of J^^^^ . 

We now proceed to show that jC^^^ approaches uniformly to J^^*^ as e goes to 0. Since the infinity 
is an attractive fixed point the Julia set must be in a finite region of the complex plane. We assume that 
they are inside of the disc of radius R, i.e., \z\ < R, z E J*^^*^. Therefore we can bound |-E(2)| by 



\E{z)\ < sVr^I zeJ^^^'. (52) 
The summation of ( ^l| ) can be estimated as 

2 \p'E (A^-'-'BX) I < ( £ \pA 3VrVW\ = ^-^3Vi?^, (53) 



fc=0 \fe=0 
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which vanishes as |e| approaches to for all integer s because we assume — jA'!" < 1. This proves 
that all points in J^^^^ approach to J^^^^ uniformly in the integrable limit. 

In the above we considered the case of |A| > 1, |A'| < 1. The other case |A| > 1, |A'| < 1 can be also 
treated similarly if z; is transformed into wi = 1/z; in (p3|). Since this transformation is equivalent to the 
exchange of the role of A and A' (and replacement of 6 hy —6 ) in (p3|), we can replay on the ty-plane the 
same argument to the above. 

We conclude this paper by showing pictures which represent the convergence of the Julia set to the 
points of iterative maps of the integrable system. The parameters of the map are fixed at A = 4 and 
9 — O.OStt. Under the choice of these parameters, A', hence e, can be changed freely. In the integrable 
limit e = 0, J^^*^ = {-4^e"* "3'^", n = 0, 1, 2, • • ■}.lt is shown in Fig. 3c. As e differs from the Julia 
set expand from these points as seen in other pictures. The real and imaginary axes are not drawn except 
in Fig. 3a, so that the points in the neighbourhoods of z = —4 and are visible in other figures. 



Fig. 3a 



Fig. 3b 
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Fig. 3c 
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By studying analytical property of a piece of Toda lattice we attempted to clariiy how a nonintegrable 
system approaches to the integrable one. Our argument is based on the fact that the two dimensional 
Toda lattice can be disjoined into small pieces, which are integrable by themselves and are called Toda 
molecules. A Toda molecule is composed from smaller pieces, which we called Toda atoms. Hence the 
two dimensional Toda lattice is a crystal consisting of Toda atoms. For such a macroscopic system being 
integrable every piece must be joined very carefully not to create a Julia set. 

In this onnection it will be worth while recalling that a similar property is possessed commonly in 
other integrable models. In the solvable lattice models the partition function is factorizable into a product 
of Boltzmann weights. The Yang-Baxter equation is a condition imposed on the factors to be connected 
properly. Another exampele is the factorizability condition imposed on the string amplitudes which led 
us to the r function of the KP hierarchy. In any case the connection rule must be such that the symmetry 
characterizing the unit blocks is preserved under the coupling. 

We will be interested in studying analytically properties of the compound system of two GLM pieces 
in the forth comming paper. 
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